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Abstract 

Linear theory of the parametric beam instability or the self-amplification of parametric x-ray radiation 
(PXR) from relativistic electrons in a crystal is considered taking into account finite emittance of the 
electron beam and absorption of the radiation. It is shown that these factors change essentially the 
estimation of threshold parameters of the electron bunches for the coherent X-ray generation. The 
boundary conditions for the linear theory of the effect is analyzed in details and it is shown that the 
grazing incidence diffraction geometry is optimal for the growth of instability. Numerical estimations of 
amplification and coherent photon yield in dependence on the electron current density are presented for 
the case of mm-thickness Si crystal and 100 MeV electrons. Possible improvements of the experimental 
scheme for optimization of the coherent radiation intensity are discussed. 


1 Introduction 

The advent of X-ray free electron lasers (XFEL) opened new era in the investigation of matter on Angstrom 
lengthscale with fs time resolution muHii]. However, need for GeV electron accelerators and hundreds 
meter long undulator modules results in high construction and maintenance costs. At the moment there 
are two operating facilities urn which are extremely overbooked and several projects under construction 
1318]. In recent years there were investigated several mechanisms that could lead to compact lab-size bright 
and coherent x-ray sources ini nnuni US- In this contribution we will theoretically analyze the possibility 
to achieve the x-ray lasing from 100 MeV electrons in mm-thick crystals based on the parametric beam 
instability effect. This effect was predicted by Baryshevsky and Feranchuk in 1983 m, they showed that 
above threshold current density value the interaction between parametric x-ray radiation electric field and 
relativistic electrons leads to instability and exponential growth of radiated intensity. This effect was realized 
in the THz range with artificial periodic structures la, but for x-rays and crystals the threshold current 
density was estimated to be 10^A/cm^. Such current density values were considered to be unrealistic at 
the time of publication m, but became recently available from short electron bunches. It is also important 
to stress that such bunches could pass through the mm-thickness crystal without its destruction during 
the production of x-ray radiation m, the concept being close to ’’diffraction before destruction” one, used 
intensively in the XFEL imaging approaches [16]. 

Before quantitative estimations, let us qualitatively describe what we expect of beam instability. In case 
of incoherent PXR the intensity of radiation is proportional to the number of electrons in the beam, as they 
are arranged randomly in the beam. However, emitted field affects trajectories of the electrons, this process 
may be called self-interaction. It may happen that rearrangement of the electrons takes place and they 
become ordered in space due to self-interaction. In that case electrons start to radiate coherently and one 
can expect radiation intensity to grow up to proportional to squared number of electrons in the beam. Eor a 
nC bunch this may result in 10^^ increase of intensity compared to spontaneous PXR case, the origin of this 
enhancement is the same as SASE process in XEEL. The aim of the present contribution is to investigate the 
parameters of the electron beam and experimental geometry under which the coherent amplification starts. 


1 


entry surface 



exit surface 


Figure 1: Coordinate system conventions. 


In the present contribution we use the classical electrodynamics and linearized theory that enable to 
describe the initail phase of intensity amplification and locate the necessary for the parametric beam insta¬ 
bility parameter area. We take into account absorption in the crystal, consider growth of instability out of 
incoherent PXR under grazing incidence geometry and estimate influence of the electron beam’s emittance. 
The paper is organized as follows. The theoretical basics is developed in Sec|^ The development is done 
in following steps: In Sec |2.2| we calculate the the electric current originating from an external filed. It can 
be represented as JE {in k — w domain), where J is a 3 x 3 matrix,it is the operator J that describes 
self-interaction. In Sec |2.3| we derive the dispersion equation in case of beam instability in order to be able to 

we analyze the boundary conditions and interplay between spontaneous and amplified fields. Based on the 
analysis, the amplitudes of out-coming waves are calculated. Sec {^provides several numerical examples. In 
Secj^the obtained results are analyzed and ways for improvement are outlined. 


find the field eigenwaves in the crystal. Sec 12^ describes how to pick out the most amplified waves. In Sec 2.5 


2 Calculations 

2.1 Notations 

In present paper we make use of Gaussian-cgs units. Further, we try to avoid frequency w itself and use 
the combination w/c instead. Hereafter we will imply w/chy w in order to shorten equations. The beam 
electrons’ velocity is described by dimensionless value j3 = v/c. The coordinate system we use is shown if 

Figj^ 

n is a surface normal, |n| = 1, gf is a vector of reciprocal lattice space, it is parallel to the crystal’s surface 
as we consider grazing incidence geometry. In this problem we will have to deal with two strong waves k 
and kg = k ^ g. k will either lie in xz-plane or be close to it. 

2.2 Homogeneous part of the electron beam current 

In this subsection we consider two forms of the beam current density, they correspond to zero and nonzero 
emittance of the electron beam, in case of nonzero emittance we regard distribution over directions of 
electrons’ velocities in the beam only. Although distribution over magnitudes of velocities has influence on 
the effect, impact of distribution over directions is more significant. Let us first consider the zero emittance 
case. The beam is assumed to be uniform. We represent coordinate and velocity of an electron of the beam 
numbered j as 


= j3 + dj3j{t) 

Vj {t) = Vjo + (3ct + 5r{t) (1) 

Srj{t) = cS(3j{t) 
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( 2 ) 


Sl3j{t) is obtained from the differential equation which is consequence of the second Newton’s law pT] : 

df3j = ——{{1 - f3j <S) Pj)E(rj,t) + Pj X H{rj,t)), 


here 7 ^ = 1 /y 1 — iSj^SjSj and Vj are functions of time, m is an electron’s mass, ( 8 ) denotes the outer 

(Kronecker) product. In the framework of linear approximation S/3j is found from Eq.Q under substitution 
/3j —> /3 in the right-hand side. Then, the homogeneous part of the beam current can be found by the 
following expression 


j{r,t) = -ecY^Sl3j{t)S (r - rj{t)) (3) 

j 

In this approximation after routine calculations one obtains 


ie^Ub , / 3 ( 
j{k,w) = - ( 1 -' 




mc^w ' w — k • jS {w — k • 

here - is the concentration of electrons in the beam ( 1 /m^), 


/3 ® (5)E{k,w) 


(4) 


j{k,w) = y j{r^t)eyi^i{wt — k-r)drdt (5) 

The approach to calculation of the formula @ is described in m pp. 146-153] and [18] in detail. Also, 
Eq.Q can be derived by a relativistic approach like in [191 pp. 144-145]. One can obtain an expression for 
current in a frame of reference, where mean velocity of electrons of the beam equals zero, and then transfrom 
into the laboratory frame of reference by using Lorentz transformations. 

In present research we consider the vicinity of Cherenkov’s resonance, that means \w — k • j3\ « w. It 
lets us omit first and second addends in Eq.Q. Eurthermore, for the wave kg we can neglect the whole 
beam current by the same reasoning. 

As to nonzero emittance case, then we make use of the approximate expression for the current and take 
into account distribution over directions of velocities by performing an integration 


j{k,w) 


le 


mc^w 


/ 


driB 

df3 {w — k • 


fd 0 f3df3E{k^ re). 


( 6 ) 


here |/3| in the integration area is implied to be constant. Before we choose the integration region we have 
to introduce a vector ko, which is an initial approximation to the solution of dispersion equation so that k 
can be represented as k = ko n5z^ where 1^;^ | « w. In this case w — k' j3 = {w — ko • (5) — ^ z- The form 
of the velocity distribution does not change the qualitative behavior of the function (|^, and we will choose 
the from which makes the equations easier. It was convenient for us to choose a rectangular integration 
area in parameters {w — ko • /3z. We denote X = —{w — ko • f5) for compactness. Also, we assume that 

dnB/{dXdf3z) = const in the integration region. Visualization of such integration area is given in Eigj^ . 
We chose this particular distribution because it can be integrated analytically. However, it is realizable as it 
does not have any singularities or other peculiarities except for its abrupt border. Still, one can consider it 
as a model of real distribution with realistically sharp border. 

Although exact expression for the integral in Eq.([^ was found, it is quite cumbersome and there is no 
need to present it here. We are much more interested in asymptotic behavior of this integral at Sz 0. 
Concerning this, it was established that at arbitrary limits of integration over X the integral tends to a 
constant when Sz 0 , but if either X^in or Xmax equals zero, the integral acts like {l/Sz)B when Sz 0 , 
where B is some operator independent on Sz- Erom now on we consider the case X^in = 0? even though 
results are essentially the same in either case. Here is an expression for B 


B = 


ie^UB k‘^ — 

TTiC^W Xrfjidx 


Bn 

0 

Bi3\ 


0 

B 22 

0 

(7) 

Bsi 

0 

B 33 J 
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Figure 2: Integration area in coordinates Pz- lies in PxPz plane. In more details, the area is placed 
on the surface of a sphere |/3| = const and the figure represents a projection on the plane PxPz- 


where 


— ^2 ^^Ozi + ^kQz^Pzmin Pzmax)) W ^ 

B22 — j ^2 2^oiPzmin Pzmax) ('^ ^ 

-^33 7\^Pzmax “ 1 “ Pzmin) -^31 -^13 ~i 7\^^z^Pzmin Pzmax^^ 

Z Kqx ^ 

"’■‘J 

koy = 0, for coordinate system see Figj^ Fortunately, it turns out that under certain conditions the 
operator B can be represented in a form of a diad proportional to /3 0 /3 with a relative error ^ 10“^, where 
/3 is one from such /3 that w — ko • ^ = 0. We’ll consider such configurations only. Finally, coming back to 
universal notation, one can state that in nonzero emittance case 


In 13 z 


-Pz 


- ln/?z 


2.3 Dispersion equation 

In the framework of two-wave approximation Maxwell’s equations lead to the following linear system(see [T71 
pp.146-153]), where we have already omitted negligible parts of the current; 

(T - fe (g) k)E - w'^X-gEg = -Y(3 ® (3E 

{Tg - kg<S) kg)Eg - = 0, 

where 


wi n o 7— i a\2 zero emittance case 

Y = (10 

nonzero emittance case 

here = 47re^n5/(mc^), C is a coefficient dependent on parameters of distribution of electrons’ velocities 
over directions, kg = E = E{k, re). Eg = E{kg^w)^ eo is a mean susceptibility of the crystal, x^, X-g 

are Fourier components of the crystal’s susceptibility, T = k‘^ — w‘^eo^ Tg = kg — w‘^eo. 
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In general, there is a straightforward way to obtain the dispersion equation, it is to find the determinant 
of the 6x6 matrix corresponding to the system of equations ^ and equate it to zero. One may even 
take more accurate expression for current. However, here we present more demonstrative derivation of 
approximate dispersion equation for cr-polarization in case of zero emittance which proves to give solutions 
almost identical to those of the straightforward dispersion equation. 

Let us regard such /3 that can be presented as 

/3=^fe + /3x 

with appropriate accuracy. Where /3|| = • k) /k^ f5±_ \\ k X kg \\ Bfj - the cr-polarization unit vector. 

In this case 

(i.E= ■ k + 

k 

Scalar products of Eqs. © and Ccr give 


TE- - w\_gE- = (^E • k + 

TgEg - uP'XgE'^ = 0 

We don’t know E • k yet. Scalar products of Eqs.(|^ and k, kg give 

{T-k^){E.k)-w\.g{Eg.k)=-{(3-k)Y(^^{E.k) + p^E<^) 

T{E-kg)-{k- kg) (E-k)- w^x-g {Eg ■kg)=-{l3- kg) Y {E ■ k) + ^^E^^ 

"^9 {^9 ' k^ ~ {k • kg^ (^Eg • kg^ “ w‘^X q (-^ • k) = 0 

(Tg - kl) {Eg . kg) - w\g {E ■ k) = Q 

This is a linear system of equations for the variables {E • k) ^{E • kg) , [Eg • k) , [Eg • kg). Thus, one 
finds [E • k) through and can substitute it into Eq. ©• The latter becomes a homogeneous linear system 
for and Eg. In order that it has non-trivial solutions the determinant of the corresponding matrix should 
be equal to zero, which yields the following dispersion equation for cr-polarized waves 


{w-k-(3)^D^ = ^ _ ;32 ^ 


(13) 


where D^j = TTg—XgX-g"^^ - left-hand side of XRD(X-ray diffraction) cr-polarization dispersion equation, 
here the zero-emittance case was considered. Eq.([T^ is valid in the vicinity of intersection of Cherenkov 
roots and one from diffraction roots. The addend D^j on the right-hand side of Eq.(13) is due to non-zero 
E • k. It was verified that its role is significant. Hence, the longitudinal component of the field cannot be 
neglected and the approximation of transverse waves fails. 


Eq.(13) can help to get insight into the situation, to find a range of better parameters etc. In order to 
obtain an exact dispersion equation corresponding to the linear system § we used the following technique. 
Substituting Eg from the first equation of Eqs.(|^ into the second one we get 


[Tg - kg® kg)[[T - k ® k) T YjS ® (5) - w^XgX-g 


E = {) 


(14) 


Scalar products of this equation with vectors k^kg^fS results in a homogeneous linear system for the variables 
{E • k) ^ [E • kg) ^ [E • (5) . If one equates the determinant of this linear system to zero, then one obtains the 
dispersion equation corresponding to Eqs. <§• Erom our viewpoint this approach is preferable, because it 
does not make us deal with any coordinate systems. Resultant expression contains scalar products of the 
vectors k^kg^fB only. Obtained dispersion equation is presented below. Y is defined in Eq.(lO). 
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Figure 3: Intersection of dispersion surfaces. Green sphere - \k\ = w, brown sphere - \kg \ = w, blue plane - 
w — k • ^ = 0 - the Cherenkov plane. It is supposed that /3 = const and w = const here. Thick purple line 
denotes the points where Bragg’s condition is fulfilled. 


= Y 


D^{TgW^€o{G^ - w^eo/S^) - W{M^ - Tg/S^)} - V%W 


where G = k • 13, M = kg • 13, V = k • {kg X /3), W = XgX-gW^^ 


(15) 


~ (^0 XgX-g) 


TTg - W 


/ k‘^k‘1 - {k ■ kgf\ 

\ w\el-XgX-g)l 


- left-hand side of the XRD dispersion equation for 7r-polarized waves. 

If in Eq.(15) right-hand side equaled 0, it would be the XRD dispersion equation for total electric field. 
It would yield 8 diffraction roots. In considered case of beam instability {w — k • f5)‘^ or {w — k • f5) in 
denominator of Y increases order of the equation by 2 or 1 respectively. The diffraction roots remain almost 
unchanged, additional 2 or 1 roots are by few orders of magnitude less than the diffraction roots in our 
parametrization, we will call them Cherenkov roots. Here by roots we mean Sz in k = ko Y n5z- We will 
number all roots of Eq.(15) by m = 1..10 or 1..9. m = 1..8 correspond to diffraction roots, the rest 
- Cherenkov roots. Until this moment we distinguished the two forms of current by corresponding beam 
emittance. Now we see, that we might as well call them two- and one Cherenkov roots cases. It will be used 
further in Secl2.5[ 


2.4 Intersection of diffraction and Cherenkov roots 

Our aim is to maximize the increment of parametric beam instability. Hence, we have to choose an appro¬ 
priate region in k — w for investigation. It corresponds to intersection of Cherenkov and one from diffraction 
roots. Eirst approximation to solution of this problem is shown in Eigj^ 

The surfaces in EigJ^are situated in a vector space of k^kg^g. The two spheres are the same as in the 
problem of two wave dynamical diffraction theory (see [20l p.l31]). The purple line corresponds to fulfillment 
of Bragg’s condition(see [20l p.l21]) and noticeable interactions of the waves k and kg. Apart from this 
intersection in case of bunch instability there should be an intersection of this line with the Cherenkov plane. 
To be precise, there is no exact intersection with the Cherenkov plane in Eigj^ we can only speak of the 
point of the plane closest to the purple line. However, as already said, it is only first approximation, we will 
deal with the vicinity of this ’’intersection” in more detail below. 
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Figure 4: Relative orientation of ko, f3 in case of intersection of one from diffraction roots and 
Cherenkov roots. The angle 4> is defined by d. The blue circle represents the freedom in choice of (3. 


The region of parameters under investigation can be formulated in other words: Bragg’s condition (2k • 
= 0) should be approximately fulfilled and the angle between k and /3 should be small. 

Need for intersection of roots is seen in Eq.(13). The term (w — k • fd) and , consequently, the increment 
is larger as the term D^j is smaller. 

Let us discuss how to find more exact intersection and first approximation to solution of the dispersion 
equation ko (the solution is /c = /cq + We have to consider ko with real components only in order to 

satisfy intersection with the Cherenkov plane w — ko • [3 = ^. We introduce an angle 0 between ko and (3 
and rewrite last equation as 


w = koP cos6>, 


( 16 ) 


from now on let us consider ko and 0 as free parameters, then Eq.(16) sets the value of w 


Eurther, in order to achieve intersection with one from diffraction roots, ko should be a root of either Dcr 
or Dj.. We chose the former case. That is 


(fcg - w^eo){klg - W^€o) - XgX-gw"^ = 0, 

where w is calculated by Eq.(16). One can represent kog in the following way 


(17) 


^ 0 ^ 


kl ^2ko • g ^ g‘^ = ko ^ kid, 


(18) 


where d describes deviation from Bragg’s condition. 


Owing to Eqs.(16) and (18) Eq.(17) turns into 


(1 — eo cos^ 6>)(1 + d — eo cos^ 0) — XgX-g cos^ 0 = 0, 


(19) 


which is a linear equation for d. However, exact d obeying Eq.(19) is complex as eo, XgX-g complex. 
It contradicts the statement that the components of ko are real. There may be several ways to overcome 
this obstacle, we replaced eo, XgX-g Eq.(19) with their real parts and solved it for d. Then, with good 
accuracy we can say that we are as close to intersection of diffraction and Cherenkov roots as possible. 

Having set values of ko and 0, calculated by the above described procedure parameter d sets relative 
orientation of ko and g, whereas the angle 0 corresponds to relative orientation of /3 and ko, although there 
remains certain freedom in choice of orientation of /3 (it is shown in EigQ. There is similar freedom in choice 
of ko with respect to g as well. 

Even though in the above paragraphs we tried to find solution of the problem step by step, now we can 
say that we merely have a system of equations relating g,j3,ko and w corresponding to the intersection. 
Hence, if one sets certain relative orientation of g,l3 and n in experiment, then one can approximately know 
the most amplified out-coming waves. 

It is necessary to note that if imaginary parts of eo, XgX-g vanished, then there would be exact intersection 
of one from diffraction roots and Cherenkov roots. In the vicinity of intersection D^j would be proportional 
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to {w — k • (3) or, in other words, to dz- Left-hand side of Eq.(13) or (15) would be proportional to 5^ 
(in zero emittance case). When calculating Cherenkov 6 z from Eq.(13) or (15) it is enough to retain the 
term and substitute (5;^ 0 in all other positions as Cherenkov 6 z is several orders of magnitude less than 

diffraction roots. Thereby, in case of zero imaginary parts of the crystal susceptibilities (see [m pp.146-153]) 
Cherenkov 6 z ^ ^ where j is the electron beam’s current density. Clearly, the same behavior holds 

for increment. In present research we retain imaginary parts of cq, XgX-g thereby take into account 
absorption. In this case in the vicinity of ’’intersection” is not proportional to 5z^ it is approximately 
constant, although the constant being quite small. That fact reduces order of 5z in down to in zero 
emittance case and down to 5z in nonzero emittance case. To sum up, under absorption in zero emittance 
case dz ^ in nonzero emittance case dz ^ j. 


2.5 Boundary problem 

We use the following approach to estimation of the growth of instability. All but one electrons of the beam are 
considered as a medium with certain susceptibility, whereas the remaining electron is assumed to travel along 
the crystal unperturbed by the field, providing inhomogeneous part of the current. Eor more information on 
this method see I211122112I] and 121 pp.377-402]. The equations for the field become 


(T — /c (g) k)E — w^x-g^g = ^ Hf36 {w — k • f3) 

(Tg -kg0 kg)Eg - = 0, 


( 20 ) 


here H = Sir'^iwefc, the addend with the Dirac delta function is the unperturbed electron’s current, see for 
example [22]. One can obtain equations for the field in case of incoherent PXR by substitution T ^ 0 in 

Eqs.([^. 

General solution of Eqs.(20) is the sum of general solution of the homogeneous system of Eqs. and a 
particular solution of the inhomogeneous system of Eqs.([^. In can be expressed like (see [21]) 

( 21 ) 


E{k, w) = E\k, w) + Y. (^« - 


here i and h stand for in- and homogeneous respectively. E'^{k^w) is derived from Eqs. (20) in assumption 
that 5 — k^\w)^ is an ordinary variable, kn - normal to the surface component of the wave vector, 

k^T'^ - m-th solution of the dispersion equation for kn- The directions of E^{w) are determined by Eqs.(|^ 
at corresponding ki^\ while their magnitudes are arbitrary constants. One may apply similar consideration 
to the field in vacuum. Then the standard boundary conditions are (see m) 


/ 


Ec{k^ + kl,w)dkl = / E,{k^ + kl,w)dkl 


( 22 ) 


here c and v stand for crystal and vacuum respectively, by E here we imply total field, the sum of in- and 
homogeneous contributions. In general, there may be any other continuous across the boundary function in 
place of E{k^w) in Eq.(22). It was proven that E^ its normal derivative, the electron beam charge density 
and current comply with Eq.(22) with a good accuracy. 


In present contribution we consider grazing incidence geometry (see [20l pp. 154-160]), such choice will 
be justified below. Under this geometry all the diffraction roots should be taken into account. In the case 
of diffraction only, there are 8 /c-waves and 8 kg-wdives. By k- or kg-wdives here we mean waves with in¬ 
surface components of wave vector equal to those of ko or ko+9 respectively. Amplitude of each kg-wdive is 
unambiguously related to amplitude of one from A;-waves. Hence, one can speak of only 8 eigenwaves of the 
field in the crystal, half of the waves being direct waves and the other half being specularly reflected ones. 

As to the problem with beam instability, as in the previous case, the amplitudes of kg- and A;-waves 
are related for both diffraction and Cherenkov roots. Moreover, the relation remains the same, because the 
second equation of Eqs.([^, relating the amplitudes, is identical to that in case of X-ray diffraction. Hence, 
the unknowns are 10 or 9 (depending on the form of the beam current) scalar amplitudes of the eigenwaves 
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Figure 5: Sketch of the out-coming waves. T, Gu^ Gd are amplitudes of the waves shown by green and 
brown arrows, k^- and kg^ are in-surface components of wave vectors of the out-coming waves. 


and the vector amplitudes R^ T, Gu^ Gd (actually, they have only two independent components as the field in 
vacuum is transverse) (see Figjs]). Summing up, one has 18 or 17 unknowns respectively. The straightforward 
way to formulate the boundary condition problem is to write down the equations of continuity of the field 
itself and its normal derivative on both boundaries of the crystal plate, moreover, one has to include into 
the system of equations the continuity of the beam current and charge density on the entry surface. More 
precisely, to get a well-conditioned linear system one should take the projections of the equations for the 
field and its derivative on corresponding polarization vectors from similar problem for incoherent PXR. It 
produces 16 equations. Further, in case of one Cherenkov root these equations should be supplemented by 
continuity of the beam current (although current density is a vector, the three equations are proportional 
to each other and only one independent equation is left). In case of two Cherenkov roots apart from the 
current continuity equation one can make use of continuity of the beam charge density. Hence, one has 18 
and 17 independent equations for the cases of two and one Cherenkov roots respectively, which is enough 
to resolve corresponding boundary condition problems. Such consideration was performed and is exact in 
the framework of present research. Later it was noticed that there is much more demonstrative approximate 
solution to the formulated boundary problem. It undergoes only small deviations from the exact solution at 
sufficiently low beam currents. Fortunately, such currents coincide with real electron beam currents. It is 
this approximate solution that is to be described below. 

First, we’d like to state one from results of this approximate consideration. Let us deal with one Cherenkov 
root case and thick crystal plate(/ >> LsraggExt.^ here LsraggExt. is the secondary extinction length). 
Suppose one has a solution for the boundary condition problem for incoherent PXR. Now, if one wants to 
get -R, T, Gu^ Gd in the case with bunch instability, then he only has to multiply T and Gd by ex.p{i6i^h) 
and leave R and Gu unchanged, is the Cherenkov root. Thus, it is obvious why there is an exponential 
increase in intensity of radiation from this perspective. 

Let us prove it. Both in- and homogeneous contributions to the field are expressed through the Dirac 
delta functions and after integration in Eq.(22) there are nonsingular amplitudes left. From now on we will 
denote by E these integrated over amplitudes and will formulate the equations for them. 

It is interesting that if one derives E'^ from Eqs.(20) then one finds that the inhomogeneous field itself 
vanishes. It was proven analytically for both variants of current. There were obtained expressions f /2 • (re — 
k-f5)‘^5{w — k-f5)dkn and j fi'{w — k-f5)5{w—k-f5)dkn for the cases of two Cherenkov roots and one Cherenkov 
root respectively, /i ,/2 - certain functions with no singularities. However, in the equation for continuity of 
current the contribution due to the inhomogeneous field persists and equals exactly HfdjPzi which is the 
integrated over k^ current of the unperturbed by the field electron. The contributions are non-zero because 
they are calculated by expressions of the following form, j j 2 • {w — k • /{w — k • 5{w — k • l3)dkn and 

J ji ' {w — k • ^)/{w — k • f5)5{w — k • fd)dkn for the cases of two Cherenkov roots and one Cherenkov root 


9 






respectively, ji, j 2 - certain functions with no singularities. Thus, the contribution to the current due to the 
inhomogeneous field is in the origin of the whole effect. 

Let us now consider the case of one Cherenkov root only. There is Sz in the denominator of the expression 
for the beam current density. Magnitudes of the diffraction roots are much bigger than that of the Cherenkov 
root. It was verified by performing numeric calculations that the beam currents due to diffraction eigenwaves 
can be neglected. The above arguments let us finally write down the equation stating continuity of the 
electron beam current density 

JEg = Jo (23) 


,where J = —Y(3 ® /3, JEg - the contribution due to the Cherenkov eigenwave , Jq = H/3/Pz - the 
contribution due to the inhomogeneous field. Eq.(23) doesn’t contain amplitudes of diffraction eigenwaves. 
Hence, the Cherenkov eigenwave’s amplitude can be found immediately through Eq.(23). The remaining 
diffraction amplitudes and R^Gu^Gd can be found from equations of continuity of the field and its 
derivative 


entry surface 
'^Em = R 

m 

Y,k^YEm = -w,R 

m 

Em — Gji 

m 

^ ^ ^QmEm ~ 


exit surface 

Y,exp{iSi^n)Em=T 

m 

= w^T 

m 

Y,eMi5^Yl)9mEm = Gd 

m 

= w.gGd 


(24) 


,where m = 1..9, m = 9 corresponds to the Cherenkov root, the rest - diffraction roots, Wz = ^/w‘^ — 

Wzg = Qm = Xgiw^^^o “ feg(4™^) <8) feg((5i'"^))/(eoTg((5i'"^)), = koz + I is a thickness 

of the crystal slab. Again, one should take projections of Eqs.([^ on polarization vectors from the analo¬ 
gous problem of incoherent PXR. One may notice that the boundary problem for incoherent PXR can be 
formulated in a similar way, namely, by substitution Eg E'^, -> 0, - inhomogeneous field in case of 

incoherent PXR. Eurthermore, now we will see that Eg ^ E'^. Erom Eqs.(20) one has the following relation 
for Eg (it is valid for any Em at corresponding 


,which follows from Eqs. ©• 


DEq = j Eg 


w 


n = (T-k^k)-— 1- 


kg ( 8 ^ kg 
w^eo 


(25) 


(26) 


In Eq.(25) D is calculated at 6 z = Actually, one can put = 0 in and still preserve a good accuracy 
as 5^z^ is much smaller than kzo and diffraction roots (by magnitude). But, of course, we save 5z = S, 


- x(9) n 


denominator in J. 

Similarly, E'^ is the solution of equation 


DE^ = Jo 


here D is calculated at 6 z 


0. Together with Eq.(23) one has the following set of equations 

DEg 


JEg DE^ = Jn 


JEg — Jn 


(27) 


(28) 


here D is calculated at Sz = 0. In Eqs.(28) it’s clear that Eg ^ E'^. Thereby, we’ve developed the following 
algorithm for formulation of boundary condition problem for PXR with bunch instability. One has to write 
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Figure 6: Dependence of amplification with respect to incoherent PXR on thickness of the crystal slab. 
Blue plots correspond to zero emittance case, green - to nonzero emittance case (blue and green plots are 
obtained at different beam current densities). 


down equations for boundary condition problem in case of incoherent PXR and then replace E'^ in the 
equations corresponding to the exit surface with E'^ ex.p . Applying this approach it is obvious that 

in the limit 0 or(and) I 0 one obtains incoherent PXR. 

Eventually, if we consider a thick crystal plate(/ >> LBraggExt)^ then the equation system (24) splits in 


two independent systems (just like in case of incoherent PXR or X-ray diffraction). One for the entry surface, 
where one can neglect the specularly reflected waves, and the other for the exit surface, where the direct 
waves are negligible. The former system becomes identical to that in case of incoherent PXR, hence, R and 
Gu do not encounter any growth. The latter system is similar to that in case of incoherent PXR, but E'^ is 
replaced with E'^ exp now. Hence, the magnitudes of T and Gd in case of PXR with bunch instability 

are by a factor of | exp \ bigger than in case of incoherent PXR. Thus, there is an exponential growth. 

Now we can name a reason for considering grazing incidence geometry. We have just proved that es¬ 
sentially the enhancement is defined solely by | exp That is, geometry affects the growth through 

nothing but Cherenkov root In dispersion equation ( |l3| or (15) Cherenkov root is hidden in the term 
{w — k-/3) = —Pz^z- II is clear that as ^z decreases 5z increases. And it is grazing incidence geometry that 
corresponds to small ^z- As one changes geometry the other parts of Eq.(13) or (15) are changed too, but 
the above effect proved to be much more significant. 

In case of two Cherenkov roots the main logic remains essentially the same. But now one has to replace 
E'^ with a certain combination exp (id^z^VjE^ + exp (i5^z^^l)Eio^ where E^ + E^o = E'^. Exact magnitudes 
of Eq and Eio are unknown yet. They can be found by virtue of the equation for continuity of the electron 
beam charge density, see [25]. Routine calculations show that there should be the following substitution 


E'^ 


exp 




r(9) 




+ exp 


r(10) 


^(9) _ ^(10) 




(29) 


3 Numerical results 

Fig|^presents numerical calculation of amplitude amplification obtained by accurately formulated boundary 
condition problem. 

The Si crystal was considered, surface normal is (0, 0,1) and is parallel to z-axis (n), (h, /c, 1) = (0,4, 0). 
The parameters describing geometrical configuration in case of zero emittance are shown in Eigj^ Their 
numerical values are cp = 1°27’44.061'', Ob = 30°0’30.193", OBg = 30°0’0.629", 2'kcIw = 0.13578522nm, 
'ip = 0°17’22.336", = 0°47’31.108", the beam current density je=o = 1-25 x 

In finite emittance case, geometrical configuration is like in zero emittance case except for the direction of 
the beam electrons’ velocities, in both cases the beam electrons’ energy is = lOOMeE. In finite emittance 
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Figure 7: Rdative orientation of the wavevectors and the parameters defining it? supplementary to as 
Figj^ Here 9b = 7r/2 - Ob, Obq = 7r/2 - Obq- 



Figure 8: Angular distribution of emitted by one electron quanta {N is a number of photons, is a solid 
angle) for incoherent PXR (a) and PXR with beam instability (b). ai and a 2 are measured in mrad. Gd- 
wave is considered. 


case we have to set borders of the distribution over directio ns. Here we used 
8.07 X 10-3 


= 0, Xynaxjy^ = 3.75 X 10 see Sec 


2.2 


-^zmin — 6.34 X 10 , Pzmax — 

The normalized emittance corresponding to 
these values is ye = 1.42 x mrad assuming focusing to 0.1/im^, the peak current was taken as / = 10kA; 
the emittance and peak current values are on the frontier of achievable values (see [26l |27l EH] ) • If one uses 
approximations similar to ^ and analogue of (13), one can find that the increment depends on ratio //e^, 
i.e. beam brightness. Thereby, to merely check existence of the effect on practice one should use an electron 
beam with calculated emittance (or less) and peak current and direct it so that its mean /3z lies within the 
interval {Pzmin, Pzmax)- Then, one can expect amplified (with respect to incoherent PXR) amplitudes of the 
waves leaving the crystal through the exit surface in directions described in the above paragraph. 

In Fig shows angular distribution of emitted quanta and comparison with that under incoherent PXR, 
the case of zero emittance is considered. The same parameters as in Figj^are used except for thickness of 
the crystal slab I which is increased to 3mm. ai is an increase in from Figj^ a 2 is a decrease in Obq^ For 
considered region in ^ 1 , 0^2 we performed an integration and obtained the following estimations for photon 
yield. 1.5 x 10~^photons/electron in case of incoherent PXR and 3.7 x 10~‘^photons/eleetron in case of PXR 
with beam instability. Obtained value for photon yield of incoherent PXR (in grazing geometry) is of the 
same order as in [29] , for details of PXR under condition of grazing incidence diffraction see also [301131]. 
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4 Discussion and conclusions 


Let us now return to Fig|^ and discuss it. Firstly, we see that the waves leaving the crystal through the 
entry surface do not undergo amplification at all. It is in agreement with predictions of the approximate 
solution of boundary condition problem. 

Secondly, we chose such parameters that amplification is of the same order in both cases on purpose. 
This way we see that the character of the shift to the regime of exponential growth is slightly different. The 
plot is sharper in zero emittance case. In nonzero emittance case there is only one Cherenkov root and one 
term ex.p{i5i^h) in equations, which provides steady and smooth increase. On the contrary, in zero emittance 
case there is a struggle between two Cherenkov roots, see Eq.(19), one from which provides increase, while 
the other provides decrease. That is why there is longer plateau, short depression and sharper increase in 
the end, where one root finally overpowers the other. However, at larger energies the depression disappears, 
but the effect as a whole slowly diminishes. 

Further, after I = 1mm amplification exponentially increases and tends to infinity, but one has to 
remember that it is a linear theory of the effect. It describes solely the startup process. In order to adequately 
describe further evolution one should apply more exact non-linear theory. It will result in saturation, when 
the intensity becomes constant and proportional to where N is a number of electrons in the beam, 
instead of N in case of incoherent PXR. That is, the electron beam starts to emit as a single charge. The 
calculated photon yield for PXR with beam instability (3.7 x photons/electron) increased by two orders 
with respect to incoherent PXR (1.5 x 10~^photons/electron). That is much less than the number of electrons 
in the beam N and corresponds to the startup of the process far from saturation. It makes the estimation 
for photon yield reliable. 

As to obtained angular distributions of emitted photons, then in Figj^ one can recognize a well-known 
pair of incoherent PXR peaks (see [13 P-31]), although in this case they are asymmetric due to peculiar 
geometry (grazing incidence). In FigJ^ we see that the peak, corresponding to higher angles to the crystal’s 
surface, undergoes stronger amplification. 

In present contribution we allowed for absorption in the crystal and showed that it considerably reduces 
the increment of instability {—lm{6z)). In zero emittance case Sz ^ {j is the beam current density), in 
nonzero emittance case 6 z ^ j- For comparison, under vanishing absorption 5z ^ this case is studied 
in detail in im pp.146-153]. 

In nonzero emittance case we presented results for electron beam with constant emittance. However, as 
the beam passes through the crystal its properties change, the most considerable influence being the Coulomb 
scatter on the crystal atoms’ nuclei. If one takes into account corresponding increase in emittance during 
passing through the crystal, then the effect is essentially suppressed. We figure that the effect requires the 
beam to be restrained from expanding direction-wise. This may be fulfilled to a certain extent by channeling 
of the beam electrons. This way, considerable part of the beam electrons travel along the crystal within a 
definite range of directions, characterized by Lindhard angle 0^ (see [32],[33l p.32]). In course of time the 
amount of electrons in channeling regime descends. Still, if, for instance, the amount decreases as l/y/z 
(see [34l|35]), one can expect an integral gain proportional to \fl (that is, amplification ^ exp (const x \//)), 
where z is the in-crystal coordinate, I is the crystal plate’s thickness. The beam instability in combination 
with channeling of the electrons will be considered in further research. 

One more possible improvement is to use beforehand nano-modulated electron beam. In [26] a new 
method of obtaining such electron beams is described. Nano-modulation does not change the increment of 
instability, but it permits to reduce the length of the plateau in Figj^ In this case fast exponential growth of 
intensity may start from smaller crystal thicknesses, that also reduces the Coulomb scatter negative effects. 
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